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AESTRACT 

J This  paper  is  concerned  with  the  propagation  of  shocks  and  phase 
boundaries  in  elastic  bars.  We  consider  materials  for  which  the  one¬ 
dimensional  stress  response  is  piecewise  linear  and  not  monotonic.  In  the 
presence  of  an  applied  load  the  dynamical  fields  are  described  by  a  set  of 
functional  equations.  These  equations  are  treated  asymptotically  for  a  model 
problem  involving  a  load  which  approaches  a  constant  value.  The 
dynamical  fields  approach  the  solution  given  by  a  corresponding  Riemann 
problem  at  a  rate  ( t~n  where  n  <  -2  is  given  in  terms  of  the  stress 


response. 
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SIGNIFICANCE  AND  EXPLANATION 


When  the  stress-strain  law  governing  the  dynamics  of  elastic  bars  is 
nonlinear,  certain  loads  give  rise  to  travelling  shock  waves.  If  the  stress- 
strain  relation  is  not  monotonic,  then  these  shock  waves  may  also  describe 
phase  boundaries.  In  the  event  the  applied  load  approaches  a  constant  value, 
the  fields  in  a  semi-infinite  bar  evolve  in  time  toward  a  simple  wave  solution 
of  an  appropriate  Riemann  problem.  The  nature  of  the  approach  to  this  simple 
wave  depends  on  the  detailed  loading  program.  For  materials  in  which  the 
stress-strain  law  is  piecewise  linear,  the  solution  is  governed  by  different 
linear  wave  equations  in  distinct  regions  of  an  (x,t)-plane.  The  curves 
separating  these  regions  obey  a  set  of  functional-differential  equations. 

These  equations  are  derived  and  treated  asymptotically  in  order  to  determine 
the  large-time  dependence  of  the  solution  on  the  complete  loading  history. 

The  rate  of  approach  to  the  simple  wave  solution  is  found  to  be  independent  of 
the  early  loading  history. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


ON  THE  ASYMPOTIC  ANALYSIS  OF  TRAVELLING  SHOCKS 
AND  PHASE  BOUNDARIES  IN  ELASTIC  BARS 


Thomas  J.  Pence 


I.  INTRODUCTION.  This  investigation  is  concerned  with  the  propagation 
of  shocks  and  phase  boundaries  in  elastic  solids.  Attention  is  restricted  to 
one-dimensional  motions;  for  simplicity,  imagine  a  bar  in  which  transverse 
displacements  are  absent.  The  problem  is  introduced  in  section  II,  where  it 
is  reviewed  how  change  of  phase  phenomena  can  be  modelled  by  means  of  a 
nonmonotonic  stress-strain  law.  These  laws  have  been  studied  previously  in 
[1],  [2],  [3];  a  relevant  experimental  study  is  [4].  This  section  also 
treats  the  simple  wave  that  develops  whenever  a  nonzero  load  is  suddenly 

applied  to  the  end  of  the  bar.  This  simple  wave  would  be  expected  to  mirror 
in  some  fashion  the  ultimate  state  of  affairs  whenever  the  bar  is  gradually 
loaded  at  one  end  to  the  level  0^  provided  waves  are  not  subsequently 
reflected  back  from  the  opposite  end.  Issues  involved  in  such  an  asymptotic 
study  are  discussed  in  the  third  section.  Section  IV  addresses  special 
considerations  for  materials  in  which  the  stress-strain  law  is  piecewise 
linear.  Then,  in  section  V,  we  carry  out  an  asymptotic  analysis  for  an 
example  problem  involving  such  a  material. 


II.  FORMULATION  OF  THE  PROBLEM  AND  THE  RIEMANN  SOLUTION  FOR  IMPULSIVE 
LOADING.  Consider  a  homogeneous,  semi-infinite  elastic  bar  which  occupies 
x  >  0  in  a  reference  configuration.  Pure  longitudinal  motion  is  governed  by 
the  momentum  equation 


(2.1) 


32u  3o 


where  u  =  u(x,t)  is  the  longitudinal  displacement  of  the  bar  and  a  is  the 
stress  along  the  axis  of  the  bar.  The  density  in  the  reference  state  is  taken 
to  be  one.  Let 


(2.2) 
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denote  respectively  the  strain  and  velocity  in  the  bar.  For  elastic 
materials,  the  stress  at  time  t  at  a  location  which  was  originally  at 
position  x  is  completely  determined  by  the  value  of  e(x,t)  by  means  of  the 
constitutive  relation  a  =  0(e)  .  The  sound  speed  of  the  material  can  be 
identified  by  expressing  (2.1)  in  characteristic  form  and  is  found  to  be 

/o' ( e )  .  Here  the  '  symbol  is  the  usual  shorthand  notation  for 
derivative.  We  shall  focus  attention  on  a  hypothetical  material  for  which  the 
stress-strain  law  o(e)  is  given  by  the  smooth  curve  in  Fig.  1.  Here 


0(0)  ■  0  and  we  shall  restrict  attention  to  positive  values  of  o  and  c. 
This  entails  no  loss. in  generality  since  compressive  motions  can  be  treated  by 
Considering  a  corresponding  extensional  problem  for  a  material  with  a  stress- 
strain  curve  found  by  reflecting  the  original  curve  through  the  origin. 

Ohlike  the  state  of  affairs  in  gas  dynamics,  this  technique  gives  correct 
results  even  in  the  presence  of  shocks  [5]. 


a\c) 


Fig.  1.  Stress-strain  curve 


0(C) 


Fig.  2.  Restricted  stress-strain  curve 
indicating  the  jump  associated 
with  a  change  of  phase 


the  descending  portion  of  the  curve  in  rig.  1 • is  its  most  conspicuous 
feature.  In  an  equilibrium  setting,  strains  associated  with  this  portion  are 
found  to  be  unstable  (1].  In  the  noneqpi librium  case  the  sound  speed  is 
imaginary  ar.d  the  equation  of  motion  is  elliptic  at  these  strains.  These 
difficulties  can  be  overcome  by  precluding  these  strains.  This  is 
accomplished  in  a  natural  fashion  by  considering  an  inverse  to  the  o(e) 
relation,  namely  strain  as  a  function  of  stress.  Since  the  o(e)  relation  is 
not  monotonic,  there  are  innumerable  such  inverse  functions.  We  take  the 
particular  inverse  that  leads  to  the  clipped  or  restricted  curve  as  shown  in 
FI9.  2.  The  separated  strain  intervals  are  now  associated  with  different 
material  phases.  Here  the  value  of  the  transition  stress  is  assumed  to 

be  a  property  of  the  material  (see  (4]). Several  other  methods  for  dealing 
with  non-monotonic  constitutive  laws  have  also  been  studied  (see  (6],  [7], 
(•]). 


Zf  a  load  0g(t)  >0  is  applied  to  the  end  of  the  bar  beginning  at 
time  t  ■  0,  (2.1)  is  to  be  solved  subject  to 


(2.3)  o(c(0,t))  -  oQ(t)  ,  t  >  0  , 

(2.4)  C(x,0)  ■  0  1  v(x,0)  "0,  x  >  0  • 
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The  condition  (2.4)  indicates  that  the  bar  is  taken  to  be  initially 
undefortned  and  at  rest.  The  condition  (2.3)  can  also  be  written 


(2.5) 
where 

curve  in 'Fig.  2,  thus 


e(0,t)  =  eQ (t)  , 


t  >  0  , 


where  s^ft)  is  the  strain  associated  with  the  stress  o^tt)  by  means  of  the 


(2.6) 


o(eft(t))  =  on(t)  . 


In  the  event 


a  At) 
0 


exceeds  the  value 


V 


e0(t) 


will  be  discontinuous 


and  the  strain  field  e(x,t)  will  necessarily  include  a  discontinuity  front 
associated  with  a  change  of  phase.  In  addition,  other  shock  discontinuities 
of  a  more  familiar  kind  may  arise  from  the  intersection  of  characteristics 
associated  with  different  sounds  speeds  of  the  nonlinear  0(e)  relation. 
Across  any  such  discontinuity  front,  say  x  =  s(t),  the  jump  in  field 
quantities  are  to  be  restricted  by  the  shock  conditions 


(2.7) 


H  M  +  M  -  o  ,  If  M  ♦  M  -  o 


In  general  the  problem  given  by  (2.1)  -  (2.7)  does  not  admit  simple 
analytical  solutions. 

An  important  practical  problem  associated  with  this  system  is  that  of 
impulsive  loading.  By  this  is  meant  the  situation  where  oQ(t)  is  9*ven  by 

oQ(t)  5  am.  In  this  case  the  problem,  which  is  given  on  the  quadrant 

x  >  0,  t  >  0,  is  a  variant  of  what  is  known  as  the  Riemann  Problem  [2].  The 
solution  follows  from  the  absence  of  both  length  and  time  scales  in  the 
initial  and  boundary  conditions.  It  is  given  by 


(2.8) 


e(x,t)  =  e(X)  ,  v(x,t)  =  v(X) 


where  X  =  x/t  ;  e(X)  can  be  found  by  a  construction  which  is  outlined  in 
the  following  paragraph.  Other  treatments  of  similar  problems  may  be  found 
in  [2]. 

Let  o(e)  be  the  upper  convex  envelope  on  the  interval  [ 0 , ]  of  the 
clipped  curve  o(e)  of  Fig.  2.  Here  is  the  root  of  <*(£„,)  ■  (>„.  The 

A 

curve  o(e)  will  in  general  consist  of  a  number  of  line  segments  connecting 

A 

portions  of  the  clipped  curve.  The  curve  o(e)  inherits  the 
differentiability  of  the  original  restricted  curve  o(e)  and  so  will  be 
smooth  at  all  values  of  strain  in  the  interval  0  <  e  <  en  with  the  possible 

A 

exception  of  e  *  a.  At  the  value  e  =  a,  the  derivative  c'(e)  may  -  or 

A 

may  not  -  be  discontinuous  (see  Fig.  3.)  In  the  event  o'(e)  is 

A 

discontinuous  at  e  =  a,  o’ (a)  is  to  be  regarded  as  given  by  the  interval 

A  A  A 

[o'(a-t),  o' (a-)].  In  this  way  the  graph  of  o'(e)  is  a  monotonically 
decreasing  curve  on  0  <  t  <  c  .  Hence  the  equation 


(2.9) 
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Upper  convex  envelopes  and  corresponding  simple  waves  for  two  values 
of  0^.  Here  O' (e)  is  discontinuous  at  e  ■  a  in  the  example  on  top; 
no  such  discontinuity  in  8'(C)  occurs  in  the  bottom  example. 


has  solutions  £  for  each  value  of  X  on  the  interval  [/o'U^),  ro  '  (0  )  J  . 

A 

Whenever  0(e)  contains  line  segments,  (2.9)  will  have  an  interval  of 

solutions  in  £  for  certain  values  of  X,  say  X  =  X4  ...  X  .  At  all  other 

l  n 

values  of  X,  (2.9)  uniquely  determines  a  function  e(X).  Next,  e(X)  is 


extended  to  0  <  X  <  »,  X  *  X^  by  defining  e(X)  =  0  for  X  >  /o'(0)  =  X{ 


and  e(X) 


for  X  <  /o' (£  ) . 


This  construction  orders  the  sound  speeds  of  the  rays  x  =  Xt  in  such  a 
manner  that  values  of  strain  associated  with  higher  sound  speeds  are  located 
further  down  the  bar  compared  to  those  strain  values  with  lower  sound 
speeds.  In  the  process  it  naturally  positions  shock  and  phase  boundary  curves 
x  =  X^t,  as  well  as  ensuring  the  correct  values  of  £  on  the  x  and  t 

axes.  By  virtue  of  (2.1),  (2.7),  the  function  v(X)  is  found  from  e(X) 
through 

0  X  >  X,. 


(2.10) 


v(  X )  = 


;<v  - / 


e(X-) 


'a' { s )  ds 


xi+i  <  X  <  xA  , 


where 


Xn+1  5  °> 


v(Xj) 


v(Xf)  +  XjefAj-)  -  £(Xj)]  . 


The  solution  (2.8)  of  the  impulsive  load  problem  thus  consists  of  a 
partitioning  of  the  quadrant  x  >  0,  t  >  0  into  sectors  by  rays  through  the 
origin.  Across  these  rays  the  solution  may  be  either  smooth  or 
discontinuous.  In  the  sectors,  the  strain  and  velocity  fields  are  either 
constant  or  arise  from  solutions  of  (2.9).  In  both  cases  it  can  be  shown 
that  at  least  one  of  the  two  Riemann  invariants  associated  with  (2.1) 
remains  constant.  This  being  the  case,  the  solution  is  said  to  be  a  simple 


It  is  worth  mentioning  that  if  at  some  time  t^  >  0  the  bar  is 

subsequently  impulsively  unloaded  back  to  c*0,  the  wave  pattern  which 
results  (before  any  interactions  with  the  loading  waves  which  at  t^  are 
further  down  the  bar)  can  be  found  by  an  analogous  construction  involving 
lower  convex  envelopes  of  0(e)  . 


III.  NONIMPULSIVE  LOADING.  Whenever  the  load  0Q(t)  is  not  simply  a 

positive  constant,  the  solution  of  the  system  (2.1)  -  (2.7)  is  not  a  simple 
wave  similarity  solution  like  that  of  the  previous  section.  Instead  one  must 
take  account  of  both  families  of  characteristics  associated  with  (2.1). 
Suppose,  however,  that  the  applied  lead  eventually  attains  or  merely 
approaches  the  final  value  0B  .  In  this  event,  one  expects  the  corresponding 

simple  wave  solution  with  oQ(t)  =  oa  to  give  the  ultimate  number  of  shocks 


and  phase  boundaries,  and  to  yield  the  correct  limiting  order  of  the  waves  in 
the  bar  as  t  tends  to  infinity.  The  spacing  of  these  waves  will  depend  on 
the  manner  in  which  Oq( t)  <3^,  nevertheless  the  value  of  the  dynamical 

fields  in  the  limit  A  =x/t  fixed,  t  ♦  •  is  given  by  the  simple  wave 

solution. 


The  approach  to  the  asymptotic  state  may  be  studied  by  decomposing  the 
displacement  u  into  the  corresponding  simple  wave  solution,  and  a  correction 

A 

to  be  denoted  by  u  .  By  virtue  of  (2.8),  (2.2)  the  former  is  expressed 

A 

ufl(A)t;  the  latter  u  is  assumed  to  be  o(t)  as  t  *  •»,  A  fixed. 

Thus 

A 

(3.1)  u  *  uQ(X)t  +  u(A,t)  , 

while  the  governing  equation  (2.1)  becomes 


2A  A2  2A  ^  ./111 

(3.2)  —  u,  +  —  u, . - u,  +  u  -  c»*(—  u,  J  — —  u,  %  *  0  . 

t2  A  t2  AA  t  At  tt  vt  A''  fc2  AA 

Here  subscripts  of  A  and  t  denote  partial  differentiation.  Entering 

(3.2)  with  (3.1)  one  obtains  to  leading  and  second  order 

0  -  {A  [x2  -  »•<«;>]»;}. 

+  |a_  [xj  -  -fi;u  ♦  ;tt  -  ^  i2x  -  om«;)]  } 


(3.3) 


+  terms  higher  than  second  order. 


In  arriving  at  the  above  result  it  is  necessary  to  make  use  of  the  expansion 

ux)  -  o’fu’)  +  o-(uj)  ux  . 

The  expression  in  the  last  parenthesis  of  (3.3)  dominates  all  succeeding 
expressions  for  large  times  and  so  must  independently  vanish.  Thus  one  draws 


(3.4)  ufl(A)  “  +  c2 

for  constants  c^  ,  C2  »  or 

(3.5)  Uq(A)  »  J*  a’”1(s2)ds 

where  the  superscript  denotes  functional  inverse.  It  is  easily  verified 

that  (3.4)  gives  rise  to  the  constant  strain  and  velocity  fields  of  the 
impulsive  load  problem,  while  (3.5)  corresponds  to  solutions  of  (2.9).  The 
Riemann  solution  discussed  in  the  previous  section  is  the  unique  way  in  which 
it  is  possible  to  assemble  these  solutions  (3.4)  and  (3.5)  in  a  manner 
which  satisfies  (2.4),  (2.7),  o0<t)  =  om  and  also  jumps  over  the  specified 

Interval  (a,0). 


-6- 


The  nature  of  the  approach  to  this  simple  wave  is  governed  by  the 
expression  in  the  second  parenthesis  of  (3.3).  Requiring  this  expression  to 
vanish  leads  to 


(3.6) 


-21 


UAt  +  "  Utt  =  ° 


in  the  case  of  (3.5),  while  for  (3.4)  it  is  easiest  to  revert  back  to 
independent  variables  x  and  t  whereupon  one  finds  that  the  equation  can  be 


written 


u. .  =  o' (c,)u 
tt  1  xx 


The  latter  has  the  familiar  solution 


(3.7) 


u  =  a(x  -  /o' (C^)  t]  +  b(x  +  /o'(c^)  t) 


for  any  functions  A( • )  and  B( • ) .  The  former  (3.7)  is  also  easily  solved 
(3.8) 


a  1 

u  =  “  C(/X  t)  +  5 ( A )  =  t  C(xt)  +  D ( x/t ) 
/A 


for  any  functions  C ( * )  and  D(«).  Unfortunately  the  boundary  condition 
(2.3)  cannot  be  used  directly  for  determining  the  free  functions  which  have 
emerged  from  this  treatment.  Instead  one  anticipates  appropriate  conditions 
to  arise  from  a  more  penetrating  analysis  of  various  short-  and  intermediate  - 
time  solutions  to  the  problem,  each  of  which  is  appropriate  in  a  different 
region  of  the  (x,t)  -  quadrant.  In  general  one  expects  to  continue  the  fields 
between  these  as  yet  unknown  regions  through  a  detailed  matching  layer 
analysis.  Rather  than  pursuing  this  program,  a  method  appropriate  to  certain 
special  materials  will  be  introduced. 


IV.  PIECEWISE  LINEAR  STRESS  RESPONSE  We  turn  now  to  a  class  of  model 
materials  in  which  the  stress-strain  curve  consists  of  a  number  of  linear 


segments.  Such  models  are  often  associated  with  the  theory  of  plasticity, 
nevertheless  here  we  shall  assume  that  loading  and  unloading  follow  the  same 
stress-strain  curve. 


The  fields  arising  from  an  impulsive  load  can  be  found  by  the  procedure 
discussed  in  section  2.  As  before,  each  line  segment  in  the  upper  convex 

A 

envelope  o(e)  of  0(e)  is  associated  with  a  front  across  which  the  strain 

A 

suffers  a  discontinuity.  For  the  materials  now  under  consideration,  o(e) 
will  consist  of  nothing  but  line  segments.  It  is  convenient  to  distinguish 
among  three  types  of  discontinuity  fronts  by  drawing  a  distinction  between  the 

A 

line  segments  comprising  o(e).  We  shall  say 

a  phase  boundary  is  a  front  which  is  associated  with  a  line  segment 

which  spans  a  clipped  portion  of  the  original  curve 
a(e), 

a  contact  discontinuity  is  a  front  which  is  associated  with  a  line  segment 

which  coincides  with  the  original  curve  o(e),  and 


a  shock 


is  a  front  which  is  associated  with  a  line  segment 
which  neither  coincides  with  the  original  curve,  nor 
spans  a  clipped  portion  of  the  curve. 


The  solution  of  the  impulsive  load  problem  will  consist  of  partitioning  of 
the  (x,t) -quadrant  into  sectors  of  constant  strain  and  velocity.  These 

sectors  are  separated  by  either  shocks,  phase  boundaries  or  contact 
discontinuities . 

We  now  inquire  as  to  what  is  the  relation  between  this  simple  wave 
solution  and  the  fields  which  would  be  produced  if  the  bar  were  gradually 
loaded  to  the  level  0^7  As  before,  the  simple  wave  solution  will  yield  the 
limiting  order  of  waves  in  the  bar  as  t  tends  to  infinity.  It  will  also 
produce  the  correct  final  number  and  ordering  of  shocks  and  phase 
boundaries.  The  contact  discontinuities,  however,  are  not  associated  with 
true  curves  of  discontinuity  in  the  field  variables.  Instead  they  indicate 
nondispersive  wave  packets  across  which  the  dynamical  fields  gradually  change. 

The  approach  to  this  simple  wave  can  be  examined  by  exploiting  the 
linearity  of  the  material  in  the  different  strain  intervals.  The  difficulty 
lies  in  determining  the  regions  in  the  (x,t) -plane  in  which  the  strain  takes 
values  in  the  individual  intervals.  The  boundaries  of  these  curves  must  be 
either  curves  of  constant  strain  or  curves  across  which  the  strain  jumps 
between  values  from  different  intervals.  In  either  case  these  conditions  lead 
to  functional  equations  when  the  strains  are  expressed  in  terms  of 
D'Alembert's  solution  to  the  linear  wave  equation. 

In  order  to  illustrate  this,  we  consider  the  material  of  Fig.  5.  The 
particular  form  of  this  stress-strain  curve  is  motivated  by  its  similarity  to 
the  material  previously  introduced  in  Fig*  1.  We  shall  suppose  that 

lim  oQ(t)  ”  am  where  is  as  indicated  in  Fig.  5.  The  corresponding 

t-*" 

impulsive  load  solution  is  also  given  in  Fig.  5.  This  solution  consists  of  a 
phase  boundary  and  a  contact  discontinuity  separating  regions  in  which  the 
dynamical  fields  are  constant.  Suppose  further  that 

oQ(0)  -  0  ,  oQ»(t)  >  0  . 

We  shall  let  x  =  s(t)  denote  the  phase  boundary,  while  x  -  q(t)  shall 
denote  the  curve  upon  which  the  strain  has  value  |i,  thus 

(4.1)  e(q(t),t)  ■  v  . 

Let  Aq ,  A<| ,  a2  denote  the  regions  in  the  (x,t) -plane  in  which  the 
strain  lies  in  the  respective  intervals  (0,o),  (fl,y),  (y»«»).  In  each 

^2^  2  ^2^ 

region  hi  the  strain  obeys  the  equation  — —  “  c  — —  which  implies  that 

3t2  1  3x2 

(4.2)  e(x,t)  -  f^x  +  c^t)  +  gA(-x  +  Cjt) 

for  as  yet  undetermined  functions  f  and  g.  The  velocity  in  must  then 
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be  v(x,t)  =  c^fj^x  +c^t)  -  c^g^(-x  +  c^t).  The  fj  and  g^,  along  with 

s(t)  and  q(t),  are  eight  unknown  functions,  each  of  a  single  argument. 

They  must  satisfy  eight  conditions  given  by:  the  two  initial  conditions 
(2.4),  the  boundary  condition  (2.3),  the  two  shock  conditions  (2.7),  two 
conditions  stemming  from  (4.1),  and  a  condition  that  the  velocity  is 
continuous  across  x  *  q(t).  The  last  three  conditions  in  the  above  list  g' /e 
rise  to  functional  equations  since  q(t)  appears  as  an  argument  of  f^,  f2, 

g1#  g2.  Moreover,  since  s(t)  appears  in  (2.7),  the  condition  on 

x  *  s(t)  furnishes  a  pair  of  functional  differential  equations.  In  the 
expressions  which  follow,  we  shall  employ  parenthesis  solely  to  indicate  the 
argument  of  a  function. 


From  the  initial  conditions  (2.4)  it  can  be  shown  that  the  Riemann 

J  y 

invariant  associated  with  the  characteristics  -rr-  m  -c.  is  identically  zero 

at  u 

on  Ag.  Thus  f g ( z)  S  0  which,  by  virtue  of  (4.2),  (2.5),  yields 


g0U)  - 


*0(*/c0) 

0 


z  >  0  , 
z  <  0  . 


The  other  six  unknown  functions  do  not  admit  simple  solution 
representations.  In  what  follows  we  shall  restrict  attention  to  large 
times.  The  analysis  will  be  shortened  considerably  by  taking  advantage  of  the 
simple  wave  solution  depicted  in  Fig.  5.  This  is  not  necessary;  the 
asymptotic  simple  wave  could  be  deduced  from  the  analysis.  Such  a  program, 
however,  requires  the  consideration  -  too  lengthy  to  be  included  here  -  of 
numerous  possible  cases. 

V.  LARGE-TIME  DYNAMICAL  FIELDS.  The  boundary  condition  (2.5)  can  be 
incorporated  directly  into  the  functions  f 2,  92*  pi*ce  °f  the  functions 

*1*  9i*  *2'  ?2  we  sha11  instead  employ  f,  g,  h  through  the  expressions 


e(x,t)  - 

f(x  + 

c.t)  + 

g(-x  +  c^t) 

in 

A1 

(5.1) 

h(x  + 

v- 

h(-x  +  c2t)  +  eQ(t  -  x/c2> 

in 

A2 

• 

f.  - 

v(x,t)  - 

c^ftx 

+  C1 1) 

-  c^ g ( “X  +  c^t) 

in 

A1 

»v 

(5.2) 

c2h(x 

+  c2t) 

+  c2h(-x  +  c2t)  -  c2eQ(t  -  x/c2) 

in 

A2 

% 


The  five  unknowns  f,  g,  h,  s,  q  are  to  be  detrmined  from:  two  conditions 
holding  on  x  =  q(t)  which  stem  from  (4.1),  a  condition  expressing  the 
continuity  of  velocity  on  x  *  q(t),  and  the  two  shock  conditions  (2.7) 
which  hold  on  x  *  s(t).  The  first  two  of  the  conditions  holding  on  x  - 
q(t)  become 

(5.3)  f(q(t)  +  ctt)  +  g(-q(t)  +  cit)  *  y  , 


K\ 


hr. 


K* 


fc 


With  the  aid  of  these  two  equations,  the  continuity  of  velocity  condition  may 
be  written  as 


c  c 

(5.5)  +  c^t)  =  ~  h(q(t).+  c2t)  +  j  (l  -  —  Jy  . 


1  1 

The  two  shock  conditions  can  be  manipulated  into  the  form 

D1 

1 

D1 

'1 


D1 

(5.6)  [c^  +  s(t)jf(s(t)  +  c^t)  +  Yc~  ~  ®  t  sufficiently  large. 


D1 

(5.7)  [ c1  -  s(t)]g(-s(t)  +  c^t)  +  2c~  =  ®  sufficiently  large 


The  phrase  "t  sufficiently  large"  indicates  that  these  are  the  equations 
which  hold  once  the  phase  boundary  x  =  s(t)  has  become  the  leading 
disturbance.  This  eventuality  follows  from  the  simple  wave  solution,  moreover 

the  simple  wave  solution  also  indicates  that  s(t)  /o ( y ) /y  t  and 


q(t)  ♦  c 


Assume  that  the  wave  and  front  speeds  obey  the  ordering 


<5. A) 


c  >  s ( t )  >  c0  >  q(t)  >  0 


for  t  sufficiently  large.  One  may  conclude  directly  from  (5.3)  -  (5.7) 
that 


f  ( z)  *  fm  ,  g(z)  ♦  g^  ,  h(z)  ha 


as  z 


\{t)  *  1  ,  q(t)  ♦  8 


as  t  ♦ 


ft’f  values  for  fm,  gm,  and  a  are  obtained  from  a  consideration  of  the 
equations  (5.3),  (5.6),  (5.7)  for  large  times.  As  t  +  •  these  equations 


become  respectively 

f-  ♦  1- 
which  give 


D1  D1 
Cci  +  a3f-  +  it;  =  0  '  [ci  -a]*»  +  2r, 


0  , 


It  n  \  A  A  D1  yOUl) 

(5.9)  a  -  /Cj  +  —  »  / — —  , 


-D. 


-D. 


f"  ~  2c  „  [c.  +  a] 


r  g. 


2c1tc1  -  a] 


The  values 


(5.10) 


•  2  '■c  ■*'  2  c 

2  2  *  2 


’2  ' 
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follow  from  (5.4),  (5.5).  Since  q(t)  ♦  c2t  and  s(t)  ♦  * — “  t  as  t  ♦  “*, 

it  is  immediate  that  A0,  A1#  A2  tend  toward  the  sectors  of  the  simple  wave 

solution  displayed  in  Fig.  5.  Upon  entering  the  values  of  f^,  g^,  h^ 

into  (5.1),  (5.2),  the  corresponding  sector  values  for  strain  and  velocity 
are  recovered. 

In  order  to  study  the  approach  to  these  values  let 


(5.11) 


s(t)  =  t  +  Sj(t) 

q(t)  =  c2t  +  q1 ( t ) 
f(z)  -  tm  +  f  ^  ( z ) 
g(z)  =  g_  +  g.(z) 


h(z)  =  h  +  h4 (z) 

CD  1 


for  new  functions  s.,,  q1 ,  f.,,  g^ ,  h^.  Upon  substituting  from  these 

expressions  into  (5.3)  -  (5.7)  and  considering  a  balance  of  the  second  order 
terms,  one  obtains  asymptotic  expressions  for  these  functions.  This  procedure 
applied  to  (5.6)  yields 

D. 

0  -  [c,  +  «  +  ^(t)]  {f,  +  ^([c,  +  o]t  +  s.ft)}}  +  — 

»  f^s^t)  +  [c1  +  a  +  s1(t)]f1([c1  +  a] t  +  s^t)) 

~  Sl(t)fw  +  [c^  +  ajf^tc^  +  a]t)  , 


which  in  turn  implies 


(5.12) 


8,(0 


[c  +  a] 

- j - f^lc,  +  ajt) 


Similarly  (5.7)  leads  to 


(5.13) 


.  tcl  -  “I 

s^t)  -  — - - gr^C  Cc,  -  “It)  » 


while  (5.3)  leads  to 


(5.14) 


f^tc,  +  c2]t)~  -  g1(tc1  -  c2)t) 


•  -  *  .  *■  *  '•*  ■  *  ■  *  •  1 m  m  9  4. "  •*  »  "  *"  •  " 


.*-  v-  - :«,%  :•  *.  .v 


r, 


Eliminating  g.,  and  between  (5.12)  (5.13)  (5.14)  reveals  that  f1 

obeys 


'I 


r  °i +  c2~" 

.v 

1 

O 

+ 

V 

0 

1 _ 

\ 

1 

O 
— » 

+ 

O 
— 1 

(N 

O 

1 

r- 

O 

C1  ‘  ° 

’( 

°1  -  a 

l—  — J 

/ 

—  — 

\ 

L  J 

as  z  ♦ 


This  condition  will  be  written 
(5.15) 
where 


f1(k2z)  ~  k^f 1 (z) 


as  z  ♦  • 


(5.16) 


k.  = 


k1 


C1  ‘  a 


c1  +  a 


2  ^  [Cl  +  a][Cl  -  c2] 

'  k2  *  Ic,  -  a][c^  +  c2J  * 


k2  >  1  >  k1  >  0  . 


It  follows  from  (5.8)  that 

(5.17) 

The  asymptotic  expression  (5.15)  indicates  that 

(5.18)  f  ^  ( z )  ~  Azn  ,  n  «  In  k^ln  k2  , 


where  A  is  undetermined.  In  light  of  (5.17),  it  follows  that  n  <  0  so 
that  f-j(z)  is  indeed  dominated  by  f^  whenever  z  +  ■»,  as  was  assumed  in 


the  development  leading  to  (5.15).  Moreover  the  following  argument 
demonstrates  that 


(5.19) 


n  <  -2  . 


Proof:  (5.18)  implies  (k^k^)1*  =  k^1  +  ^ n,  which,  in  conjunction  with 

1  +  ^2 n  >  i .  This  last 


h 


0  <  k^k2  = 


c  —  c 

•  —  — — —  <1  and  n  <  0  yields  k. 
c  +  c  ■*  1 


1  2 

result,  along  with  (5.17),  yields  1  +  V2  n  <  0,  which  is  (5.19).  / 


One  finds  from  (5.18),  (5.12),  (5.13)  that 

-(1  +|)  •  -'n+' 


(5.20) 


2'.  n  ....  1°,  ♦«]' 

Az  ,  S,j(t)~ - jj -  At 


On  account  of  (5.20)  and  (5.19)  one  has 


(5.21) 

where  sq  is  undetermined. 


K  + «] 


n+1 


*t(tl  *  *0  -  ( n+ 1 ) f 


At 


n+1 
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The  asymptotic  behavior  of  the  corrections  h^(z)  and  q^(t)  are  found 

from  the  remaining  equations  (5.4),  (5.5).  Entering  (5.5)  with  (5.11), 
(5.10)  one  obtains 

g 

“h1(2c2t  +  q^t))  -  f1([c1  +  c2]t  +  q1  ( t ) )  =  0 
from  which  it  follows  that 

C1  C1  +  C2  n 

(5.22)  h  (z)  -  [-  ■2-c—  ]  Az"  . 

C2  2 

Equation  (5.4)  is  the  most  delicate  in  the  analysis.  The  exact 
equation  yields 

(5.23)  u  =  h^  +  h1[2c2t  +  qi  (t)  -  hl-q^t))  +  (t )/c2)  . 

In  this  equation  h1(2c2t  +  q.,(t))  =  o(1)  as  t  ♦  •.  It  is,  however,  not 
appropriate  to  expand  h(-q^(t)}  by  means  of  (5.11)  unless  q1(t)  ♦ 

Assume  for  the  moment  that  q^(t)  ♦  -*•;  then  C^-q^ (t)/c2)  ♦  e^, 
h^-q^t))  ♦  h „  so  that  (5.23)  gives  up  *=  jj.  Since  £„  >  u  the 
assumption  Is  false.  Instead  (5*23)  demands  that  q^(t)  ♦  qQ  >  with 
p  *  h^  -  h(-qQ)  +  e0(-q0/c2).  In  the  event  that  h(z)  admits  a  Taylor 
expansion  at  z  *  -qg,  then  (5.22),  (5.23)  yields 

{h(-qQ)  +  1-  e0(-q0/c2)}  [-q(t)  +  qQ]  -  h^t)  -  [c/Cj]  [c,  +  C^At"  . 

2 

This  in  turn  implies  that  q(t)  ~  qQ  +  0(tn)  whenever  the  coefficient  in 
parenthesis  is  not  zero.  Collecting  the  results  from  this  section  we  have 

f(z)  -  fm  +  Azn  +  o(zn) 


f  r* 


The  constants  fw,  hs,  gs,  a,  k1#  n  are  given  in  (5.9),  (5.10),  (5.16), 
(5.18).  The  constants  A,  qg,  Sq  are  as  yet  unspecified.  Indeed,  in  order 
to  assign  values  to  A,  qg,  Sg  one  may  surmise  that  it  is  necessary  to  take 
account  of  the  complete  loading  history  crQ(t).  Consider,  for  example,  the 
constant  Sg.  The  phase  boundary  x  *  s(t)  asymptotically  approaches  the 
line  x  =  at  +  sQ.  This  line  issues  from  the  t  axis  at  time  t  =  -sQ/a,  so 
that  Sg  defines  a  time  scale  for  the  problem.  The  only  source  of  such  a 
time  scale  lies  in  the  applied  load  aQ(t). 

The  dynamical  fields  are  obtained  by  substituting  from  (5.24)  into 
(5.1),  (5.2).  This  yields 


e(x,t)  ~  e  + 


Ac1  [c1  +  c2 

2nc  n+1 
Z  C2 


i  n 


v(x,t)  ~  -  ay  c2[eo#  -  y]  + 


I-i- * 1 1  - -t  J 

Acltc1  +  C21 


tn  , 


„n  n+1 
2  C2 


1  -K-fl 


in  region  A2.  In  region  A-|  one  obtains 


e(x,t)  ~  y  +  A 


v(x,t)  - ay  +  Ac 


c.  +  a  n+2 

t°t +  £  1"  *  1  [-,-f] 

,  n+2 

n  c.  +  a 

+  -fel  I-t  -  f  1 


t"  . 


Although  the  value  of  A  is  not  found  from  this  analysis,  the  exponent  n 
governing  the  rate  of  approach  to  the  simple  wave  is  determined  from  both  the 
stress-strain  behavior  of  the  material  and  the  ultimate  level  of  the  applied 
load  by  means  of  (5.18). 
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